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ABSTRACT 

Context. Recent observations of brown dwarf spectroscopic variability in the infrared infer the presence of patchy cloud cover. 

Aims. This paper proposes a mechanism for producing inhomogeneous cloud coverage due to the depletion of cloud particles through 
the Coulomb explosion of dust in atmospheric plasma regions. Charged dust grains Coulomb-explode when the electrostatic stress of 
the grain exceeds its mechanical tensile stress, which results in grains below a critical radius a < being broken up. 

Methods. This work outlines the criteria required for the Coulomb explosion of dust clouds in substeiiar atmospheres, the effect on 
the dust particle size distribution function, and the resulting radiative properties of the atmospheric regions. 

Results. Our results show that for an atmospheric plasma region with an electron temperature of © = 10 eV (=» 10^ K), the critical 
grain radius varies from 10“’ to lO^'* cm, depending on the grains’ tensile strength. Higher critical radii up to 10“’ cm are attainable 
for higher electron temperatures. We find that the process produces a bimodal particle size distribution composed of stable nanoscale 
seed particles and dust particles with a > with the intervening particle sizes defining a region devoid of dust. As a result, the 
dust population is depleted, and the clouds become optically thin in the wavelength range 0.1-10 pm, with a characteristic peak that 
shifts to higher wavelengths as more sub-micrometer particles are destroyed. 

Conclusions. In an atmosphere populated with a distribution of plasma volumes, this will yield regions of contrasting radiative 
properties, thereby giving a source of inhomogeneous cloud coverage. The results presented here may also be relevant for dust in 
supernova remnants and protoplanetary disks. 
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1. Introduction 

A defining characteristic of brown dwarfs is that their atmo¬ 
spheres are cool enough that dust particles can condense and 
form large-scale cloud structures (Helling et al., 2004; Dehn, 
2007; Helling & Woitke, 2006; Helling et al., 2008a,b; Witte 
et al., 2009, 2011; Tsuji, 2002; Allard et al., 2001; Burrows 
et al., 1997; Marley et al., 2002; Morley et al., 2012). Cloud 
formation depletes the ambient gas and alters the objects’ resul¬ 
tant electromagnetic spectrum (e.g. see recent review by Helling 
& Casewell (2014)). Observationally, the presence of clouds 
has been associated with photometric variability, which is be¬ 
lieved to be caused by inhomogeneous cloud coverage (Bailer- 
Jones & Mundt, 2001; Artigau et al., 2009; Radigan et al., 2012; 
Gillon et al., 2013; Girardin et al., 2013; Buenzli et al., 2014). 
However, such variability could also be connected to magnetic 
activity (Littlefair et al., 2008; Metchev et al., 2013), such as 
asymmetrically distributed, magnetically induced spots (Scholz 
& Eisloffel, 2005). Recent observations of Luhman 16B’s in¬ 
frared spectroscopic variability have allowed a two-dimensional 
map of the brown dwarf’s surface to be constructed. It reveals 
large-scale bright and dark features supporting the idea of patchy 
cloud cover (Crossfield et al., 2014). This cloud variability can 
be explained using a spot model where thick and thin cloud re¬ 
gions at different temperatures (indicating different atmospheric 
heights) provide the required inhomogeneity (Radigan et al., 
2012; Apai et al., 2013; Burgasser et al., 2014; Morley et al., 
2014). To quantify the inhomogeneity, recent observations have 
probed the horizontal and vertical structure of the atmosphere. 


identifying pressure-dependent phase variations that underlie the 
observed variability (Buenzli et al., 2012; Biller et al., 2013). 

Inhomomgeneous cloud coverage may be significant with 
respect to the transition from L-type to T-type objects (e.g. 
see Burgasser et al. (2014)). Ultracool L and T dwarfs bridge the 
divide between cool stars and giant gas planets, and understand¬ 
ing the transition from L-type dusty atmospheres to dust-free 
T-type atmospheres is a continuing problem (e.g. Tinney et al. 
(2003); Eaherty et al. (2012); Dupuy & Liu (2012); Burgasser 
et al. (2014)). In the latter case, the atmosphere may not be cloud 
free, but rather the clouds lie below the level where the atmo¬ 
sphere becomes optically thick. Differing models have been in¬ 
voked to explain the transition, including cloud depletion due 
to changes in atmospheric chemistry (e.g. Knapp et al. (2004); 
Tsuji (2005); Burrows et al. (2006)), the formation of holes in 
the cloud layers (e.g. Ackerman & Marley (2001); Burgasser 
et al. (2002); Gelino et al. (2002); Allard et al. (2013); Crossfield 
et al. (2014); Morley et al. (2014)), and due to atmospheric 
circulation (Ereytag et al., 2010; Showman & Kaspi, 2013). 
Eurthermore, a growing body of theoretical and observational 
evidence (e.g. Hallinan et al. (2008); Route & Wolszczan (2012); 
McLean et al. (2012);Ostenet al. (2009); Antonova et al. (2013); 
Berger et al. (2009)) suggests that ionization processes can be 
significant in these dusty environments producing extensive re¬ 
gions of plasma (Helling et al., 201 lb, 2013; Bailey et al., 2014; 
Stark et al., 2013). This implies that substeiiar atmospheres 
could be composed of dusty gas-plasma mixtures, underlining 
the importance of atmospheric dusty plasma processes and their 
role in cloud evolution. 
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A plasma is defined as a quasineutral gas of charged and neu¬ 
tral particles which exhibit collective behaviour (Chen, 1984). 
Locally, the entirety of the gas does not need to be ionized for it 
to behave as a plasma. As long as the collective, long-range elec¬ 
tromagnetic infiuence of the many distant charged particles dom¬ 
inates over short-range, binary interactions the gas only needs 
to be weakly ionized (i.e. partially ionized, fe > 10“^) for it 
to exhibit plasma behaviour (Diver, 2013; Fridman, 2008). In 
substellar atmospheres, Alfven ionization can produce localized 
volumes of plasma with degrees of ionization, fe, ranging from 
10“^ to as large as 1, depending on the atomic and molecular 
species that is ionized (Stark et al., 2013). Alfven ionization oc¬ 
curs when a stream of neutral gas (that reaches a critical thresh¬ 
old speed) impinges on a low-density magnetized plasma. The 
infiowing neutrals scatter the plasma ions, leaving behind a sig¬ 
nificant charge imbalance of electrons that are accelerated by 
their collective self-repulsion to ionising energies (Alfven, 1960; 
Diver et al., 2005). In principle, the volumes of atmospheric 
plasma created can range from small to large scales, potentially 
populating a significant fraction of the atmosphere. 

Helling et al. (2011a) investigated the effect of dust-induced 
collisional ionization in substellar atmospheres using Drift- 
Phoenix model atmosphere results. In this process, if the en¬ 
ergy imparted from one dust grain to another during a col¬ 
lision exceeds the grain’s work function, an electron will be 
liberated. Turbulence-induced dust-dust collisions were found 
to be the most efficient dust-dust collisional ionization process 
that could enhance the local degree of ionization, but on its 
own it was insufficient to form a plasma. The degree of ion¬ 
ization can be further enhanced by atmospheric electrical dis¬ 
charge events (Helling et al., 2011b, 2013; Bailey et al., 2014; 
MacLachlan et al., 2009; MacLachlan et al., 2013): an ensem¬ 
ble of inter-grain discharges, in analogy to laboratory microdis¬ 
charges (Becker et al., 2006), or large-scale cloud lightning will 
greatly increase the electron number density. Furthermore, the 
bombardment of the atmosphere by cosmic rays (Rimmer & 
Helling, 2013) will also contribute to the formation of atmo¬ 
spheric plasmas. 

The properties of the produced plasma, such as the elec¬ 
tron and ion temperature, will depend on the ionization process 
which created it. Due to their significantly lower mass relative 
to the neutrals and ions, the plasma electrons only lose a small 
fraction of their energy following a collision. As a result, the 
electrons and the heavier species will not necessarily be in ther¬ 
mal equilibrium and can have differing temperatures. In typi¬ 
cal terrestrial discharges Te > Tgas ~ T,, where Ti is the ion 
temperature and electron temperatures are Te ~ 1 - 100 eV 
(=s 10"^ - 10® K) (Fridman, 2008; Diver, 2013). In substellar at¬ 
mospheres electrons from thermal ionization are in thermal equi¬ 
librium with the ambient neutral gas {Tg ~ Tgas) and have elec¬ 
tron temperatures of Tg ~ <9(0.01 - 0.1 eV) (ss <9(10^ - 10^ K)). 

Contemporary atmospheric models of brown dwarf and ex¬ 
oplanetary (substellar) atmospheres do not exploit collective 
plasma processes (such as dust growth via plasma deposition) 
and how these effects fit into a global understanding of cloud 
formation and other atmospheric phenomena (such as large- 
scale hows). A major challenge is to model cloud formation and 
its consequent effect on the surrounding environment, involv¬ 
ing plasma processes such as non-spherical dust growth (Stark 
et al., 2006) and electrical discharge events (MacLachlan et al., 
2013). Opik (1956) investigated the electrostatic disruption (or 
Coulomb explosion) of dust grains, where the electrostatic stress 
of a body holding a net charge exceeds its mechanical ten¬ 
sile strength (the maximum stress or pressure that a mate¬ 


rial can withstand while being pulled), resulting in it fractur¬ 
ing. See also Rhee (1976), Hill & Mendis (1981), Fomenkova 
(& Mendis (1992), Rosenberg & Mendis (1992), Gronkowski 
(2009), among others. 

In this paper we outline a model to contribute to the un¬ 
derstanding of inhomogeneous cloud coverage in substellar at¬ 
mospheres via the Coulomb explosion of cloud dust particles. 
Consider an atmosphere with plasma regions populated with 
dust cloud particles. The spatial extent and distribution of the 
plasma regions will depend upon the ionization process that pro¬ 
duced them and environmental factors. For example, inhibited 
plasma transport perpendicular to the ambient magnetic field can 
result in anisotropically shaped plasma regions, structured by 
the field geometry. In a dusty atmosphere, cloud particles may 
grow in the plasma environments themselves, grow in the re¬ 
gion prior to plasma generation, or have migrated into the region 
via gravitational settling, wind hows or another transport mech¬ 
anism. As the dust charges, the accumulated electrostatic stress 
on the particles may overcome their mechanical tensile strength 
resulting in them breaking up and being destroyed. Electrostatic 
disruption preferentially affects smaller grain sizes and this may 
have an effect on the particle size distribution function affecting 
the thickness and the electromagnetic properties of the cloud. 
Depending on the distribution of the dusty plasma regions on 
the substellar surface, the electrostatic disruption of dust clouds 
could yield regions of contrasting cloud properties, thereby giv¬ 
ing a source of inhomogeneous cloud coverage that could be as¬ 
sociated with the observed spectroscopic variability. 

The aim of this paper is to investigate the conditions under 
which the electrostatic disruption of cloud dust particles occur 
in brown dwarf or gas giant planet atmospheres and to assess 
the contribution of Coulomb explosions to spectroscopic vari¬ 
ability. To do so we will use an example atmosphere from the 
Drift-Phoenix model atmosphere grid characterized by log g - 
3.0, Teff = 1600 K. This paper is structured as follows: in 
Section 2 we discuss contemporary substellar atmosphere and 
cloud formation models, specifically the Drift-Phoenix model; 
in Section 3 the concept of dust in plasmas and grain charging is 
introduced. Moreover, we consider the most appropriate plasma 
sheath model required to describe the charging of the dust and 
the region surrounding the dust grain; Section 4 formulates the 
electrostatic disruption criteria of spherical dust grains in sub¬ 
stellar atmospheres using the Drift-Phoenix model atmosphere 
discussed in Section 2. It investigates the effect electrostatic 
disruption has on the dust particle size distribution, compares 
the timescale for disruption to dust growth timescales and dis¬ 
cusses the observational consequences of Coulomb explosions; 
Section 5 summarizes our findings and discusses the impact of 
the electrostatic disruption of dust on the subsequent chemistry 
that occurs in substellar atmospheres. Note that unless otherwise 
stated, the equations in this paper are in SI units. 

2. Substellar model atmospheres 

The Drift-Phoenix model atmosphere code (Helling et al., 2004; 
Dehn, 2007; Helling & Woitke, 2006; Helling et al., 2008a,b; 
Witte et al., 2009, 2011) is a ID phase-non-equilibrium model, 
that describes the formation of dust and clouds in a self- 
consistently calculated atmosphere. Drift-Phoenix model atmo¬ 
spheres incorporate radiative transfer, convective energy trans¬ 
port, chemical equilibrium, hydrostatic equilibrium (Hauschildt 
& Baron, 1999) and dust cloud formation (Woitke & Helling, 
2003, 2004). The cloud model numerically solves a system of 
differential equations that describe the time-dependent dust for- 
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mation process, where seeds form (Ti 02 ) from a highly su¬ 
persaturated gas. Once the initial seed particles have formed, 
further growth proceeds via surface chemical reactions grow¬ 
ing a grain mantle. The Drift-Phoenix model atmospheres con¬ 
sider 32 surface reactions and is restricted to seven participating 
solids: Ti02[s], Al 203 [s], Fe[s], Si02[s], MgO[s], MgSi 03 [s] 
and Mg 2 Si 04 [s]. Drift-Phoenix solves the dust moment and el¬ 
ement conservation equations, returning the parameterized par¬ 
ticle size distribution function f(a,z) of grain radius a at at¬ 
mospheric height z; and the chemical composition of the dust 
expressed in volume fractions Vj/Vtot, for all solids s involved 
where Vtot is the total dust volume. 

In this paper we are interested in the electrostatic disrup¬ 
tion of dust grains in substellar atmospheres. We consider a 
representative model atmosphere characterized by logg = 3.0, 
Tefif = 1600 K, applying results from the Drift-Phoenix model 
atmosphere simulations. We are particularly interested in the 
chemical composition of the dust particles (hence the material 
volume fractions) and their size distribution since this informs 
of their mechanical properties and the particle sizes affected. 
For this model atmosphere Fig. 1 shows the dust volume frac¬ 
tion Vj/Vtot and the mean grain size (a) as a function of atmo¬ 
spheric gas pressure pgas. The nucleation-dominated upper cloud 
deck (pgas ~ 10“" bar) predominantly contains small grains 
((a) ss 10^^ cm) composed of Ti 02 [s], whereas the lower at¬ 
mosphere (pgas ~ 10“^ bar) is mainly made up of large dust 
grains ((a) ^ 10“^ cm) composed of multiple species includ¬ 
ing Al20[s], Fe[s], Mg 2 Si 04 [s] and trace amounts of other solid 
materials. At mid-atmospheric pressures (pg^s ~ 10“^ bar) the 
mean grain size is (a) ^ 10“® cm and the grain composition is 
more mixed being composed mainly of MgSi 03 [s], MgSi 04 [s], 
Si02[s], Fe[s] andMgO[s]. 


3. Dust cloud particles in substellar plasmas 

In this section we discuss the charging and the behaviour of dust 
immersed in a substellar atmospheric plasma region and con¬ 
sider the appropriate model required to best describe the plasma 
sheath of the dust. 

3.1. Dust charging in piasmas 

Consider a dust particle submerged in an electron-ion plasma. 
We will assume that the particle currents drawn from such a 
plasma are the only contributors to the grain charge and we ig¬ 
nore other charging processes such as photoionization and sec¬ 
ondary electron emission. This assumption is valid in the ab¬ 
sence of significant UV radiation, which would liberate the elec¬ 
trons from the surface of the grain, possibly leading to positively 
charged grains (Shukla, 2001). Due to the greater mobility of 
the electrons relative to the ions, the dust becomes negatively 
charged initiating the formation of a plasma sheath (an electron 
depleted region) around the grain. As the negative charge on the 
grain builds up, the number of electrons having the appropriate 
kinetic energy to overcome the grain potential, striking its sur¬ 
face, decreases. Ions are subsequently accelerated from the bulk 
plasma, through the sheath and are ultimately deposited on the 
grain surface altering its shape, size, mass, charge and hence the 
potential of the grain. As a result the number of electrons reach¬ 
ing the surface increases further altering the surface charge. This 
occurs until the grain reaches an equilibrium state where the flux 
of electrons and ions at the grain surface is equal. At this point 
the grain reaches the floating potential (pf and the sheath length is 


of the order of the plasma Debye length T/j = (eokBTelinee^))^^^. 
The Debye length is the spatial length scale at which exposed 
charge is screened by the plasma. In this particle-flux equilib¬ 
rium state we have a negatively charged grain that is screened 
by the bulk plasma such that on length scales greater than the 
Debye length, the electric field from the grain is zero. The spa¬ 
tial extent of the sheath (the sheath length, d^h) occurs where the 
electric potential energy of the electrons is approximately equal 
to their thermal energy and the electric field from the grain is 
largely shielded from the rest of the bulk plasma.. 

For conducting grains the captured charge will distribute it¬ 
self on the surface in the lowest energy configuration such that 
the electric field inside the conductor is zero, otherwise charges 
see a net force and move seeking equilibrium. For spherical dust 
grains this leads to the uniform distribution of electric charge on 
the grain surface and the electric field can be simply described. 
However, in ultra-cool substellar atmospheres dust grains are in¬ 
sulating and so once the electrons strike the surface their trans¬ 
port is inhibited and they can not easily reach the lowest en¬ 
ergy configuration. However, in a plasma the dust grains will 
be constantly bombarded by an isotropic flux of electrons and 
ions which allows the lowest energy configuration to be obtained 
even though electron transport in the dust grain is inhibited. For 
spherical dust, the electric field from the grain will be the same 
as that from an equivalent conducting grain allowing the same 
formulae to be used to describe its electric field. 


3.2. Coiiisions in piasma sheaths 

The effect of collisions in the plasma sheath defines how the 
sheath will behave and how it should be described. In this sec¬ 
tion the significance of collisions in the sheath is quantified. The 
plasma sheath surrounding a dust particle is considered colli¬ 
sionless if the collisional length scale (the mean free path, /Imfp n) 
for ions with the neutral species is greater than the sheath ex¬ 
tent, dsh- d = t/sh/Tmfp,n < 1; Otherwise, the sheath is colli¬ 
sional (Lieberman & Lichtenberg, 2005). In other words, if the 
ions do not participate in neutral collisions while transiting the 
sheath then 6 < 1; if the ions experience multiple collisions 
then 6 > 1. In the case 5 < 1 the effect of short-range bi¬ 
nary collisions is negligible in comparison to long-range col¬ 
lective phenomena and the sheath is well described by solu¬ 
tions to the Vlasov equation (Bouchoule, 1999). If the collisional 
scale length for ions is less than the sheath length then a colli¬ 
sional model describing the sheath is preferable. Assuming that 
the sheath length is approximately equal to the Debye length, 
dsh ~ do and dmfp.n = (wnCFn)^', we can write, 


/ eoksTg 

\ rieC^ 


1/2 


( 1 ) 


where rie is the plasma electron number density, n„ is the neutral 
species number density and cr„ is the corresponding collisional 
cross-sectional area. Note that in Drift-Phoenix, Ugas refers to 
the number density of all species and n„ refers to the number 
density of the neutral species only, such that Ugas - n„ + ng + tii - 
tin + 2n(,. Thus, the degree of ionization takes the following form, 
fe - nelitin+tie) - nc/(ngas-ne)-For ion-neutral interactions, we 
assume cr„ ss nr^ where ~ 10“* cm is the radius of a neutral 
species; therefore, 




( 1 -/.) 



( 2 ) 
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where Wgas is in units of [ctn“^]. The ratio of the sheath length 
and the ion-neutral mean free path, dni, (Equation 2) is plotted 
in Figure 2 as a function of atmospheric gas pressure, pgas, for 
an example substellar atmosphere characterized by logg = 3 
and Teff = 1600 K (see Section 2). Note that the atmospheric gas 
pressure pgas, in this plot and in subsequent plots, acts as a proxy 
for atmospheric height and in general does not reflect the true 
local atmospheric gas pressure. This is a result of varying the 
degree of ionization from the values calculated in Drift-Phoenix 
due to thermal ionization, to reflect additional non-thermal ion¬ 
ization processes that occur. As a consequence the local atmo¬ 
spheric gas pressure would change and deviate from its origi¬ 
nal Drift-Phoenix calculated value. Further to this, note that for 
fe- 1, dni = 0 since n„ —> 0 and Tmfp n —> oo. For low degrees of 
ionization {fe ~ 10“^) the neutrals will be the majority species 
and the transport of the plasma ions will be inhibited by frequent 
collisions with the neutrals as they travel through the partially 
ionized medium. If the degree of ionization increases and more 
of the gas is ionized, the neutral number density decreases and 
the likelihood of an ion-neutral interaction is smaller relative to 
the previous case. In the case of a fully ionized region of plasma 
(/e = 1) there will be no neutrals for the ions to interact with and 
the sheath will be collisionless. For atmospheric gas pressures 
Pgas ~ 10“^^ - 10 bar in the model considered, the collision¬ 
less sheath criterion is satisfied, dni < 1, for all cases of fe and 
Te that characterize the plasma regions in the atmosphere. At 
higher atmospheric gas pressures (pgas ^ 10“^ bar) only plas¬ 
mas with degrees of ionization fe > 10^“* satisfy dni < 1. In 
plasma regions where 5ni > 1, the sheaths are more accurately 
described by collisional models; however, dni > 1 is only ob¬ 
tained for pgas > 10“^ bar, which is a very small fraction of 
the atmosphere; hence, we are justified assuming a collisionless 
sheath. 


3.3. OML and Bohm theory 

Following Section 3.2, if we assume a collisionless sheath and 
a dust particle of radius a, the parameter f - a/d^h defines the 
plasma sheath regime and the appropriate theory required to de¬ 
scribe it: Bohm or Orbital Motion Fimited theory (OMF) (Allen 
et al., 1957). Bohm theory (f » 1) assumes that all ions that en¬ 
ter the sheath reach the dust surface and contribute to its charge, 
whereas OMF theory (^ <K 1) takes into account that some ions 
may be deflected and fail to reach the particle surface. Assuming 
for the sheath length that c/sh ~ do we have, 


Calculating the floating potential using Bohm theory gives an 
upper limit to the actual potential value. In the regime ^ 1 

the sheath length is much greater than the particle size. In this 
scenario, it is more likely that ions entering the sheath with a 
random velocity component may not collide with the particle 
surface since its surface area is much smaller than the sheath’s. 
Instead, ions will have non-zero angular momentum and may 
participate in stable orbits around the dust particle or be ejected 
from the sheath on parabolic orbits. In this case the floating po¬ 
tential, 0/ is reduced and is calculated using OMF theory and is 
given by the expression (Bouchoule, 1999), 



Tetnijy ksTij' 


(6) 


It is important to note that the floating potential in both cases 
depends on the bulk properties of the plasma. The parameter 
^ = a/d^h (Equation 4) is plotted in Figure 3. Predominantly 
^ < 1 indicating that OMF theory is best suited to describ¬ 
ing plasma sheaths in substellar atmospheres. Since the neutral 
number density increases with atmospheric gas pressure pgas, the 
electron number density also increases with pgas for a given de¬ 
gree of ionization. As a result, the Debye length decreases with 
increasing electron number density and f increases in value as 
the atmospheric pressure increases. Holding the electron num¬ 
ber density fixed, increasing (decreasing) the electron tempera¬ 
ture has the effect of increasing (decreasing) the Debye length 
and decreasing (increasing) 

There is no universal consensus on the dynamical theory of 
ions in plasma sheaths and the resulting effect on dust charg¬ 
ing in plasmas; therefore, we present results from both Bohm 
and OMF theory in order to give a pragmatic representation. 
Bohm theory gives a valuable upper limit to the floating poten¬ 
tial whereas OMF represents an indication of how the non-zero 
angular momentum of ions in the sheath affect the charging of 
the dust particle. For the results presented here, the difference in 
the floating potentials calculated using the Bohm and OMF is 
never more than a factor ss 2. 


4. Inhomogenous cloud coverage through the 
Coulomb explosion of dust clouds 

This section calculates the criteria for charged dust grains to 
Coulomb explode and the effect this has on the resultant cloud 
particle size distribution function in substellar atmospheres. 




/ 2 \l/2 

a^] 

\eokBTej 


ss 0.1 a 




(3) 

(4) 


where in the final expression iigas is in units of [cm“^] and a is 
in units of [cm]. For ^ » 1 the sheath extent is much smaller 
relative to the particle size and so the surface area of the dust 
particle is effectively the same as the sheath surface area. In this 
regime, all ions that are accelerated from the bulk plasma and 
through the sheath reach the particle’s surface contributing to its 
charge and potential. For ^ » 1 the sheath is called a Bohm 
sheath and the floating potential is given by (Bouchoule, 1999), 


4.1. Criteria for eiectrostatic disruption of sphericai dust 
grains 

Dust immersed in a plasma becomes negatively charged through 
the collection of plasma electrons. During the charging pro¬ 
cess described in Section 3 the dust can be electrostatically dis¬ 
rupted if the resultant electrostatic stress acting on the dust grain 
(due to the captured electrons) overcomes the mechanical ten¬ 
sile strength (the maximum stress that a material can withstand 
while being pulled) of the particle, breaking it apart. Consider a 
conducting spherical dust grain of radius a with a uniform tensile 
strength F (measured in Pascals, Pa). Electrostatic disruption oc¬ 
curs when the electrostatic stress (eoF^/2) on the grain exceeds 
its tensile strength (Opik, 1956), 


(j)f 


ksTe 


■In 


Imrif 


(5) 


2 


2 \47reofl^/ 


> 2;, 


(7) 
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therefore, using Q - Ne, 

Ne > (8) 


This expression states that there is a maximum number of net 
charges that can reside on the surface of a spherical dust grain of 
radius a without it breaking up. Equivalently, this expression can 
be stated in terms of the surface charge density cr > 
which is independent of the grain size. For the critical charge 
density ctc, grains with a smaller surface area can only toler¬ 
ate a smaller number of charges in comparison to larger grains. 
Therefore, in a plasma where spherical dust particles are charged 
to the floating potential ipf (— Ea), the result can be recast in 
terms of the critical cloud particle size, below which the 
dust particles are electrostatically disrupted, 


2 \ a ) 


> s, 


(9) 


hence. 


a < a 


Coul 

crit 



( 10 ) 


Therefore, dust particles that are charged and have a size below 
a < will Coulomb explode and be broken up. For multi¬ 

composite materials (as expected for cloud particles in substellar 
atmospheres) the Rule of Mixtures (RoM) can be used to approx¬ 
imate the resultant tensile strength S of a hybrid material: 


s = m,2(y,/ytot)2„ (11) 

5'=1 


where VJViot and Sj are the volume fraction and tensile strength 
of the s* material of the composite (Callister, 1985; Torquato 
et al., 1999; Hsieh & Tuan, 2005; Kim et al., 2001; Fi et al., 
1999; Torquato, 2000; Raabe (& Hangen, 1995). The RoM can 
be used to approximate properties of composite materials such 
as the elastic modulus. Additional physical effects can be incor¬ 
porated into the RoM formulation such as interfacial bonding be¬ 
tween the composite materials; and the effect of particle size dis¬ 
tributions and particle clustering within the composite (Fi et al., 
2001). To quantify the effect of imperfect bonding between the 
individual components and the departure from perfectly formed 
idealized crystalline solids the parameter nii, e (0,1] is intro¬ 
duced. 

Using the values of tensile strengths listed in Table 1 and 
the Drift-Phoenix volume fractions y^ (Figure 1), Figure 4 
shows the critical cloud particle size, calculated using 

Equation 10 and using the tensile strength, S, from Equation 11. 
The critical cloud particulate size is plotted as a function of lo¬ 
cal gas pressure pgas for electron temperatures ranging from 1 
to 100 eV. The top, middle and bottom subplots show for 
mb = 1, 10“^ and 10^“* respectively. For mb = 10^“*, electron 
temperatures Tg - Tgas (~ 10^ K), green) are considered 

but are not considered for the cases mb = 1 and 10“^, since the 
critical particle sizes are too small and are unphysical. 

The critical particle size is dependent upon the floating po¬ 
tential and the tensile strength of the dust cloud particle, 

0// VE. For both sheath models (Bohm and OML) the floating 
potential depends upon the electron temperature, Tg. As Tg in¬ 
creases the plasma electrons become more energetic and so a 
greater number of electrons can reach the surface of the grain, 
overcoming the electrostatic potential and increasing the ulti¬ 
mate potential reached when the particle-flux equilibrium con¬ 
figuration is achieved. For the Bohm model the floating potential 


(and hence has a simple linear dependence upon Tg and so 
an increase in Tg is reflected in 0/ and 

For example when mb = 1, at pgas ~ lO^'* bar the critical 
particle radius varies from ~ 10“^ cm when Tg - \ EW 

10“* K) to 10-® cin when Tg = 10 eV 10^ K). 

Similarly the critical radius increases to ~ 10“^ cm when 
Tg - 100 eV (ss 10® K). In comparison to the Bohm model, the 
floating potential and hence calculated with OMF theory 
is smaller for a given Tg\ however, the difference in 0^°“* for the 
two models is never more than a factor ss 2. 

For a given electron temperature, the variation of with 
Pgas is predominantly determined by the variation in the total ten¬ 
sile strength of the dust particles, E, hence the volume fraction 
Ui/Utot and individual tensile strengths of the constituent mate¬ 
rials. As an example to elucidate this point, we look at the values 
of (3^™* for mb = 1 and the Tg = 10 eV (ss 10® K), Bohm theory 
solution (red-dash) in Figure 4. At the bottom of the atmosphere 
(Pgas ~ 10“^ bar), Figure 1 shows that Al 203 [s] is the domi¬ 
nant dust component {VsjViot ~ 1) which has a relatively high 
tensile strength (241.4 MPa) and this raises the overall mechan¬ 
ical tolerance of the dust. As a result, the critical particle radius 
(oc E“b 23 ^ relative low value, a®"™* ss 10“^ cm. However, as 

we climb through the cloud (pgas ~ 10^“* bar) the volume frac¬ 
tion of Al 203 [s] decreases drastically and other materials, such 
as MgSi 03 [s], that have a significantly lower tensile strength 
(9.6 MPa) dominate the mechanical properties of the grains, 
therefore lowering the electrostatic stress that can be mechani¬ 
cally tolerated. As a consequence 0^°“* increases in value from 
that at higher atmospheric pressures to ~ 10“® cm. Finally, 
at the top of the cloud the grains are predominately composed of 
Ti 02 [s], which has a tensile strength similar to magnesium sil¬ 
icates (36.41 MPa) and so the critical particle size is similar in 
value ss 10“® cm. 

cnt 

The material tensile strengths used (Table 1) are obtained 
from material science, where they are obtained under optimum 
laboratory conditions using idealized samples and present an up¬ 
per bound to their actual value. In reality the tensile strengths of 
the dust particles will be much lower reflecting the imperfect 
bonding and structure of the particles. Although the dust par¬ 
ticles could be crystalline in nature (see Helling & Rietmeijer 
(2009)) the possibility exists that they are not and are only a 
loose aggregate of weakly bonded material. The tensile strength 
is the macroscopic quantity that represents the effect of the mi¬ 
croscopic bonds holding the atoms of the solid together. An es¬ 
timation of the tensile strength of a dust grain can be made by 
considering the work done, W, by the bonds to maintain a stable 
solid structure; E ss W/Vd, where Vd - Ana^l3 is the volume of 
the dust grain. If a solid is held together by Nb bonds of energy 
Eb, then W ^ EbNb- The number of bonds can be estimated by 
calculating the number of participating atoms of radius ro in a 
cross-sectional slice through the dust grain of area hence, 
Nb ~ {alroY, where ro ~ 10“* cm is the classical radius of 
an atom. Therefore, the tensile strength of the dust grain can be 
estimated as E ss (fl/ro)^£’b/(47rfl®/3). For a crystalline spheri¬ 
cal dust grain of radius a ^ 10^“* cm and typical lattice energy 
Fb ~ 10 eV, the estimated tensile strength is E ss 10 MPa, which 
is consistent with the tensile strengths listed in Table 1. However, 
for a similarly sized dust grain that is weakly bonded, held to¬ 
gether by van der Waals bonds (Eb ~ 0.01 eV), where only 1% 
of the atoms participate in bonding, the tensile strength can be 
as low as E ss 10 MPa. Therefore, in comparison to the typical 
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tensile strengths in Table 1, E ss (9(1 - 100 MPa), the parameter 
mb can be as low as lO^”*. 

The effect of imperfect bonding is demonstrated in Figure 4, 
for values of m* = 10 “^ (middle subplot) and 10 “'^ (bottom sub¬ 
plot). Consider a^°“* calculated for Tg - 10 eV (ss 10^ K) using 
Bohm theory, at pgas ~ lO^'* bar, assuming a perfectly bonded 
crystalline dust particle {nib - 1), ss 2 x 10“® cm. For 
nib - 10 “^ this value of the critical dust radius increases by a 
factor m^*^^ ss 10 to become 0 ^°“* ~ 2 x 10“^ cm. Similarly, 
for nib = 10 ^“*, increases by a factor of ss 100 such that 
ss 2 X 10^“* cm. For aTg - 100 eV plasma the critical par¬ 
ticle radius can be as large as 0 ^°“* ~ 10 “^ cm when mb = 10 “"^. 

Charged dust grains with a radius below the critical radius 
a < 0 ^°“*, will be disrupted and broken up, leaving a popu¬ 
lation of dust grains with a > In addition, there will 

be a further population of crystalline seed particles composed 
of a single substance (e.g. Ti 02 , Fe, AI 2 O 3 , etc) with a sta¬ 
ble radius, Oseed = “ (eo/( 2 Ej))'^^ 0 y, where nib - 1 - 

For example, pure crystalline seed particles composed of Ti 02 
(E = 36.21 MPa) immersed in a plasma with Tg = 1 eV 
(ss lO'* K) will have a minimum critical radius Oseed ~ 10“^ cm. 
Consider an aggregate of crystalline seed particles forming a 
spherical particle cluster. In this context the parameter ni\, is no 
longer unity and contains a number of additional effects, such 
as the packing efficiency of arranging constituent particles into a 
desired volume. For simple cubic, body-centred cubic and face- 
centred cubic stacking geometries the packing fraction is /p ss 
0.52, 0.68 and 0.74, respectively. The packing fraction sets an 
upper limit to the value of ni\,. Assuming that the aggregate’s 
tensile strength reflects that of the crystalline seed particles, the 
ratio of the tensile strengths for the two cases is 1 /mb and the 
ratio of the critical radius is 1 / v/mb. Therefore, if the particle 
aggregate’s radius is smaller than Oseed/ ^/n%, the aggregate is 
unstable to electrostatic disruption and the unstable range is de¬ 
fined by Joseed, flseed/ 'sjnibl- However, the tensile strength of the 
aggregate may be lower than that of the constituent seed parti¬ 
cles since they may be bonded together via van der Waals in¬ 
teractions. For example, consider a particle aggregate of radius 
10 “^ cm composed of 10 “^ cm seed particles bonded together 
via van der Waals interactions (ss 0.01 eV), the tensile strength 
of the aggregate can be 10“^ MPa. In this scenario the ratio of 
tensile strengths will be smaller, resulting in a greater critical 
radius for stability and a larger region of instability. 

The seed particles that compose the particle aggregate are 
initially charged to the floating potential 4>f — 2seed/(4:^eoaseed), 
where Qseed is the charge on a seed particle. The resulting parti¬ 
cle aggregate will also be charged to the floating potential such 
that (pf - 2agg/(4:^foaagg), where gagg and Oagg is the charge 
and radius of the particle aggregate respectively. Therefore, the 
charge ratio is Qagg/Gseed = flagg/oseed- The total charge of the 
aggregate can also be calculated by summing the charge from 
the individual seed particles. The number of seed particles in the 
aggregate, Aseed, is given by the volume occupied by the seed 
particles in the aggregate, Vagg/p, divided by the volume of a 
seed particle, Vseedi therefore, Aseed = (aagg/aseed)^/p- The total 
charge of the aggregate is gtot = (aagg/aseed)^/p 2 seed, which is 
greater than the aggregate charge calculated from the floating po¬ 
tential, Qioi > Qagg- The process of particle aggregation results 
in charge being lost by the aggregate back into the plasma. 

However, this analysis assumes only surface charges; if the 
particle aggregate has a volumetric charge distribution, the elec¬ 
tric field will be non-zero inside the aggregate and the con¬ 


stituent charged seed particles will see a non-zero electrostatic 
force. This will undermine the mechanical stability of the ag¬ 
gregate and will reduce the value of nib further. This is in con¬ 
trast to a conducting dust grain, where the electric charge will 
distribute itself on the surface such that the electric field in¬ 
side is zero and the electric potential is constant. The surface 
area of the aggregate, is less than the total available sur¬ 

face area, Aseed4:^aseed- Assuming that each seed particle is iden¬ 
tically charged, the charge on the surface of the aggregate is 
given by gsurf = aseedGagg/(aagg/p) and the internal charge by 
Gint = Gagg - Gsurf- For example, assuming Oseed = 10^’ cm, 
flagg = 10 “^ cm and /p ss 0.5, the surface charge is Gsurf ~ 
lO^^Gagg- Therefore, the majority of the electric charge is held 
within the body of the aggregate. Let’s consider that the net elec¬ 
tric charge is distributed uniformly throughout the volume of 
the aggregate so that the electric field as a function of radius, 
r, inside the aggregate is £ = Qaggr j{Ane(,a\^^. Comparing the 
surface electric field and the internal electric field of the aggre¬ 
gate yields £/£surf = ^/'^agg- This allows a comparison of the 
relative strength of the internal electric field. Note that in the 
case of a conducting dust grain the internal electric field is zero. 
On length scales comparable to the seed particle size, F/Fsurf = 

aseed/Oagg ~ 10“^ for flseed = 10“^ Cm and flagg = 10“^ cm. 

Therefore, a porous dust aggregate will have a repulsive, non¬ 
zero electrostatic energy density that weakens its mechanical in¬ 
tegrity. The analysis of dust grains collected from the Jupiter- 
family comet 67P/Churyumov-Gerasimenko show the dust parti¬ 
cles are fluffy, porous aggregates with very low tensile strengths 
of the order of 10^ Pa (Schulz et al., 2015). This supports the 
tensile strength values considered in this paper and that the dust 
grains may be a loose, porous, fluffy aggregate of weakly bonded 
material. 

Therefore, in a dusty plasma region there exists an insta¬ 
bility strip for particle sizes between seed particles and parti¬ 
cle aggregates composed of multiple seed particles, defined by 
jflseed, flseed/ V^bJ- fo Order for Stable grains with a radius greater 
than flseed/ to form they cannot do so through aggregation 
nor thin film surface growth via neutral gas-phase surface depo¬ 
sition. The accretion of unstable fragments that produce a bigger 
stable particle could occur as long as the accretion timescale is 
shorter than the timescale for the fragment to Coulomb explode. 
If the plasma is only present in certain regions of the atmosphere, 
dust growth can occur classically and can be transported into the 
plasma region, whereupon the unstable particles will Coulomb 
explode upon charging to reach the floating potential. 

Electron field emission occurs when a dust grain’s surface 
electric field becomes large enough that electrons are liberated 
from its surface. If significant, the resulting current must be 
taken into account to determine the grain’s net charge from the 
steady-state balance of currents to its surface. For dust grains, 
electron field emission becomes a significant process when a 
critical electric field strength at the grain’s surface is reached 
Ffe = 4>lci ^ 10^ V cm“^ (Gomer, 1961; Axford & Mendis, 
1974; Draine & Salpeter, 1979; Ishihara, 2007; Mann et al., 
2014), where ip is the grain potential and a is the radius of the 
dust grain. This value is calculated from the Fowler-Nordheim 
equation, which determines the electron current drawn from a 
conducting dust grain with a work function, Wf (typically Wf - 
1-5 eV). Assuming (pja - 10^ V cm“*, field emission is not a 
significant grain charging process if a > \naTg — 1 eV 

plasma, (pf - 0{\ V) and so dust grains with a > 10“^ cm are 
not significantly affected by field emission, likewise fl > 10“^ cm 
and fl > 10“^ cm for plasmas with Tg - 10 eV ( 0 y = G(10 V)) 
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and 100 eV {(pf — (9(100 V)) respectively. Only for the smallest 
dust grains in a Tg = 100 eV plasma is field emission a signif¬ 
icant contributor to dust charging and so inhibit the attainment 
of for electrostatic disruption. For this grain population, 
the potential they acquire is restricted to 0 = lO^a; therefore, 
for a = 10“^ cm (10^^ cm), the potential reached is 0 = 10 V 
{(j) - \ V). As a result, the electrostatic disruption criterion for 
the grains revert to that for the case 0/ ss (9(10 V) (0/ ss (9(1 V)); 
hence, when Tg x (9(10 eV) (T^ x 0(1 eV)). Therefore, the grain 
population affected by field emission may still be able to partic¬ 
ipate in electrostatic disruption. For example, when m\, - 10 
and taking (pf - (9(10 V) (Tg - 10 eV) as the grain potential 
reached, = 10^“* cm and so dust grains with a = 10“® cm 
(i.e. those affected by field emission when Tg - 100 eV) will 
still fracture. Furthermore, in substellar atmospheres dust grains 
are primarily composed of insulating material and so the work 
function can be larger e.g. Wf — 10 eV (Miloch et ah, 2009), 
yielding larger tolerable values of E^f before field emission be¬ 
comes a significant process. 

4.2. Timescale considerations 


via hopping to a solid island of suitable composition where sur¬ 
face chemistry and chemisorption occurs, creating a new unit 
of material in the solid’s crystal structure (Helling & Woitke, 
2006). Note that the growth process is independent of the grain 
charge. Consider a seed particle in an atmospheric region that 
has just become partially ionized. If the seed particle can grow 
fast enough it may be able to surmount the critical radius 
and escape electrostatic disruption. As a result the particle size 
distribution function will not be truncated significantly by the 
process of electrostatic disruption. For a net growth velocity 
the time required for a particle to grow to a size is given by 


,Coul 

cnt 


^seeid 
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net 


(15) 


Therefore, for electrostatic disruption to have a significant effect 
on the particle size distribution function we require the charging 
time to be less than the growth time (tc < 



^r^<i 


(16) 


The electrostatic disruption of dust in substellar plasmas will 
only have a significant effect on the particle size distribution if 
the timescale for disruption dominates over other dust growth 
and destruction timescales. In a fully ionized plasma where there 
are no neutrals, nucleation and growth as described by contem¬ 
porary substellar atmosphere and cloud formation models is not 
suitable since the chemistry is predicated upon neutral gas-phase 
chemistry. In this case an extended approach of dust nucleation 
and growth in substellar plasmas is required. For a fully ion¬ 
ized plasma region, the timescale for electrostatic disruption pre¬ 
cludes the timescale for other processes associated with charged 
dust and so electrostatic disruption will occur. In a partially ion¬ 
ized plasma, the ionized regions will be composed of gas-plasma 
mixtures and so if the gas is not sufhciently ionized then the nu¬ 
cleation and growth rate for dust formation could be greater than 
the rate of destruction via disruption. 

The characteristic charging time for dust in a plasma is given 
by (Bouchoule, 1999), 
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where is the Bohm theory characteristic charging time 
for dust; is the OML theory charging time; and yo = 

-ecpfliksTe). Once seed dust particles are formed in substel¬ 
lar atmospheres they grow as they gravitationally settle un¬ 
til they ultimately evaporate at the base of the atmosphere. 
The grain-size dependent, equilibrium settling speed is obtained 
from the balance between the gravitational force and the fric¬ 
tional force exerted by the surrounding atmospheric gas on the 
dust grain (Woitke & Helling, 2003). The net growth velocity for 
heterogeneous growth is;^^"®' [cm s“'] (positive for growth; neg¬ 
ative for evaporation) which includes all the growth velocities 
for each participating condensate, = YjsXlet (Helling et ah, 
2008b). In Drift-Phoenix the growth process occurs through the 
physical adsorption of impinging gaseous molecules on the to¬ 
tal grain surface. The molecules diffuse across the grain surface 


Figure 5 shows Equation 16 evaluated as a function of atmo¬ 
spheric gas pressure using Equations 12 and 13 for the dust 
charging time and from Drift-Phoenix. Eigure 5 presents 
the worst case scenario for dust charging, where the atmospheric 
gas is just sufficiently ionized (i.e. fg - 10^^) to exhibit collec¬ 
tive plasma effects, hence Tc is at its largest; and the dust grains 
are very mechanically strong (i.e. mb = 1) resulting in hav¬ 
ing the smallest possible values and the quickest dust growth 
timescales, In general, for all cases and all atmospheric gas 
pressures < 1, so the seed dust grains charge quicker 

than the timescale for growth and the particle size distribution 
function will be significantly truncated by the electrostatic dis¬ 
ruption of dust particles. As the atmospheric gas pressure in¬ 
creases, the gas-phase neutral number density increases and so 
for a given degree of ionization, the electron number density in¬ 
creases resulting in a smaller charging time and Tg <K Eor 
higher degrees of ionization, fg > 10“^, the charging time de¬ 
creases, becoming much faster than the net growth speed of the 
dust grains and the effect of electrostatic disruption is significant 
in the evolution of the particle size distribution function. 

4.3. Effect on the dust particle size distribution function 

Cloud particles of size a < immersed in a plasma, are 

subject to fracturing by electrostatic disruption. This can lead to 
a bimodal distribution of final grain sizes, comprising particles 
larger than a^°“’ and smaller than nano-sized seed particles (both 
of which are stable to disruption), with intermediate sizes ex¬ 
punged. This could yield regions of contrasting cloud properties 
in the atmosphere, thereby giving a source of inhomogeneous 
cloud coverage. To understand the role of electrostatic disrup¬ 
tion in substellar atmospheres the effect on the dust particle size 
distribution function needs to be quantified. 

Eigures 6 and 7 show how the particle size distribution 
function f(a,z) as a function of cloud particle size a, varies 
through the atmosphere, where z is the atmospheric height. Seed 
dust particles form high in the atmosphere, the dust settles un¬ 
der the effect of gravity as it grows and then finally evapo¬ 
rates at the base of the atmosphere. As a consequence, as the 
grains grow, the peak of the particle size distribution function 
migrates towards large particle sizes as the atmospheric gas 
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pressure increases. Figures 6 and 7 shows two potential parti¬ 
cle size distribution functions derived using the dust moments, 
Kj — a^f(a,z)da, from Drift-Phoenix: a potential expo¬ 

nential size distribution (Helling et ah, 2008b) and a lognormal 
size distribution. The lognormal distribution is included to in¬ 
crease the impact and relevance of the results presented here 
to other models where such a distribution function is assumed. 
Using the Drift-Phoenix dust moments a lognormal particle size 
distribution function can be given by 


f(a,z) = -;=exp 

acr yin 


(In(a)-ju)^ 

2cr^ 


(17) 


where /r and cr are the mean and variance of the distribution re¬ 
spectively and are given by, 
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(a) - KJKo, ( 20 ) 

(A) = A 71 K 2 IK 0 , (21) 

and (A) is the mean dust surface area [cm^]. In Figure 6 and 7 
each curve represents the dust particle size distribution function 
at different atmospheric heights, starting at the cloud top where 
Pgas ~ 10 " bar and finishing at the bottom of the cloud where 
Pgas ~ 10“^ bar. The top, middle and bottom subplots show re¬ 
sults for when mb = 1 , 10 “^ and 10 “^ respectively, for the 
case Te - 10 eV (Bohm theory). As m* decreases and in¬ 
creases, larger regions of the particle size distribution f{a, z) be¬ 
come Coulomb unstable and deeper cloud layers, where grains 
are larger in size, become increasingly affected. The distribution 
function resulting from the Coulomb explosion of charged dust 
cloud particles in a plasma can be given by 
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Note that if the Coulomb explosion process manages to destroy 
all the large dust particles, all that will remain is a dust cloud 
composed of nanometre seed particles. 


4.4. Observational consequences of electrostatic disruption 

The resulting optical depth, due to the particle size distri¬ 
bution of the dust can be calculated by evaluating the following 
double integral, 

g(a,z')7Ta^ Gextdfl dz', (26) 



where Qext = 2 ext(a, 4, Vs) is the extinction efficiency of the 
dust, which is a function of the dust particle size, wavelength and 
volume fraction Vs- The optical depth is of interest since it de¬ 
termines the intensity (and therefore the flux density for a given 
solid angle) that is observed. The contrast in optical depth be¬ 
tween different atmospheric regions will be seen as a change in 
the observed flux density. Figure 8 shows the resulting dust opti¬ 
cal depth as a function of wavelength ( 0 . 1-10 pm) for a region of 
the atmosphere where the Coulomb explosion of dust particles 
occurs. To exemplify the Coulomb explosion model an atmo¬ 
spheric plasma region stretching from the bottom of the cloud 
to the top is considered, using fl^°“' for = 10 eV (Bohm) to 
calculate g(a,z) and for a lognormal distribution function. 
To calculate the extinction efficiency, a bespoke Mie scattering 
numerical routine was utilized using the Bruggeman approxi¬ 
mation (Bruggeman, 1935) to calculate the effective optical con¬ 
stants of the dust grains using values from the references listed 
in Table 2. The optical depth is calculated using a Simpson’s 
numerical integration routine over the entire dust particle size 
distribution function g(a, z) and over the entire cloud extent in z. 

As mb decreases in size and the effect of Coulomb explo¬ 
sions increase, more and more dust grains are broken up, leav¬ 
ing a thinner cloud structure with a lower optical depth with re¬ 
spect to adjacent cloud regions where the process is absent. As 
a result, the contrasting electromagnetic properties of the neigh¬ 
bouring regions would manifest itself observationally as patchy 
cloud coverage. For comparison, the optical depth for an atmo¬ 
spheric region unaffected by electrostatic disruption is shown 
(black solid line). 

For clouds with a distribution of particle sizes, dust grains 
with a length scale comparable to (or harmonics of) the wave¬ 
length of the interacting electromagnetic waves resonate most 
effectively with the radiation. The general profile of the dust 
optical depth as a function of wavelength A is determined 
by the superposition of resonant interactions between the elec¬ 
tromagnetic radiation and the dust particles. We are interested 
in photometric variability in the infrared wavelength range cor¬ 
responding to dust particles with diameters in the micrometer 
range. As a result, the optical depth peaks at a wavelength, where 
the dust cloud is optically thicker, given by Apeak ~ 2 (fl), where, 


<fl)» = 


£cZ ag(a,z)da 

^crit 

Xc™' g(a, z)da 


(27) 


Note that in the wavelength range of interest (T = 0.1 - 10 pm) 
the effect of the population of seed particles (flseed ~ 0(1 nm)) 
is negligible relative to the effect of the micrometer-sized pop¬ 
ulation, since Mie scattering (/U^ for a ^ A) dominates over 
Rayleigh scattering (A^‘^ for a <K 4). As mb decreases, the critical 
radius fl^ increases, the extent of the truncation of the particle 
size distribution increases and the resulting mean particle size, 
(fl)t, increases, shifting the peak in the optical depth to higher 
wavelengths. Further to this, since the dust particle population 
is depleted, the cloud becomes optically thinner. In the case 
mb = 10 the truncation of the particle size distribution is more 
severe with fl^°“* ~ 7 x 10“^ cm to flmax, completely reducing the 
upper cloud layers to only an ensemble of nanometre sized dust 
particles. In addition to the population of seed particles at flseed, 
the particle size distribution in lower cloud layers is restricted 
to a narrow distribution of larger dust grains. As a result of this 
narrow distribution, there is insufficient overlap in the resonant 
features at wavelengths corresponding to harmonics of the dust 
diameter, to guarantee a continuous optical depth profile; hence, 
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the dip in the optical depth at d a; 1.2 jim. This is in contrast 
to the cases mb = 10“^ - 1, where the spread in particle sizes is 
sufficiently broad that the dust-radiation resonant features form 
a continuous optical depth profile. In the extreme case where all 
dust particles fall below the critical threshold for stability, only 
a cloud composed of nanometre-sized seed dust particles will 
exist making the cloud optically thin in the wavelength range 
presented here (0.1 - 10 yum) with a characteristic Rayleigh scat¬ 
tering intensity dependence. These features present useful diag¬ 
nostics for distinguishing the Coulomb explosion process from 
other inhomogeneous cloud coverage mechanisms. 

5. Discussion 

The aim of this paper was to investigate the electrostatic dis¬ 
ruption of cloud dust particles in substellar atmospheres and to 
assess the contribution to inhomogeneous cloud coverage and 
the resulting spectroscopic variability. Dust cloud particles im¬ 
mersed in a plasma region, with a radius below a critical radius 
a < are electrostatically disrupted and broken apart. This 
results in a bimodal particle size distribution composed of sta¬ 
ble nanoscale seed particles and particles with a > with 

the intervening range of particle sizes defining a region of in¬ 
stability devoid of dust. For an atmospheric plasma region with 
an electron temperature of Te - 10 eV, the critical grain radius 
varies from s; 10"^ cm to 10“"^ cm depending on the tensile 
strength of the dust particle, characterized by mb = 1 to 10“^ re¬ 
spectively. Higher critical radii are attainable for higher electron 
temperatures. As mb decreases, the dust population is depleted 
and the clouds become optically thin. Furthermore, as the critical 
radius increases, the extent of the truncation of the particle 
size distribution increases and the resulting mean particle size 
increases, shifting the wavelength-particle size resonant peak in 
the optical depth to higher IR wavelengths. Depending on the 
distribution of the dusty plasma regions in the substellar atmo¬ 
sphere, electrostatic disruption of dust clouds could yield regions 
of contrasting cloud and optical properties, which gives a source 
of inhomogeneous cloud coverage that could be associated with 
the observed spectroscopic variability. If the entirety of the at¬ 
mosphere is ionized then all the atmospheric dust is charged 
and susceptible to electrostatic disruption. In this case inhomo¬ 
geneous cloud coverage would occur due to non-uniformity in 
the electron temperature across the atmosphere since is de¬ 
pendent upon the electron temperature. 

The electrostatic disruption of dust grains in substellar at¬ 
mospheres has a significant impact on the subsequent chemistry 
that occurs. Consider a spherical dust grain of radius a and cor¬ 
responding surface area. If such a grain is electrostatically dis¬ 
rupted it will break-up into numerous smaller fragments each 
with a smaller surface area than the initial, parent dust grain. 
However, the resultant combined surface area of the fragment 
ensemble will be greater than that of the initial dust grain. As a 
result, reagents and radicals have a greater surface area available 
to them, promoting surface catalysis. The increased collective 
surface area from fragmentation is in analogy with micronized 
coal-fired power generation: lumps of coal are broken-up into 
micrometer sized particulates (via a collisional process) the col¬ 
lective surface area of the micronized coal is greater than that of 
the parent lump. Therefore, surface catalysis is enhanced and the 
efficiency of combustion increased resulting in a greater energy 
yield in comparison to standard coal-fired stations. Although the 
specific chemical reactions in the micronized coal context are 
different, the underlying physical principle is the same: an in¬ 
creased surface area enhances the chemical reactivity. The in¬ 
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Fig. 1. Volume fraction Vj/Vtot and the mean grain size {a) of a mineral 
cloud in a substellar atmosphere with = 1600 K and logg = 3. 
Results are from the Drift-Phoenix model atmosphere simulation and 
is consistent with Helling et al. (2008b). 


creased reaction surface area may cause a stronger depletion 
of the gas-phase species that participate in grain surface chem¬ 
istry (see Table 1, Helling et al. (2008b)). Hence, localized time- 
dependent gas-phase depletion may well be caused by a sudden 
population boost of small grains induced by electrostatic disrup¬ 
tion of larger grains. 

The electrostatic disruption of dust grains may also be appli¬ 
cable in other astrophysical, dusty plasma environments, such as 
supernova remnants. Supernova remnants are associated with the 
creation of a series of shocks that propagate through the dusty 
ejecta affecting dust formation and destruction. Investigations 
have suggested that grains with a length scale < 0.05 jum are 
destroyed by sputtering in the post-shock flow while those a; 

O. 05 - 0.2 jum are trapped in the shock (Nozawa et al., 2007). 
The fraction of dust that survives the passage of the shock ranges 
between 2-20% (Bianchi & Schneider, 2007), which may be fur¬ 
ther affected by additional processes such as the electrostatic dis¬ 
ruption of dust. Furthermore, observed silicate features in proto¬ 
planetary disks are believed to be consistent with the existence 
of a population of micrometer sized grains and the absence of 
sub-micrometer sized grains (Williams & Cieza, 2011). To pro¬ 
duce such a particle size distribution, Olofsson et al. (2009) sug¬ 
gest that either dust growth processes dominate fragmentation 
at these length scales or the sub-micrometer dust grains are re¬ 
moved from the upper layers of the disk via winds or radia¬ 
tion pressure. The electrostatic disruption of dust particles could 
be an additional mechanism aiding in the removal of the sub¬ 
micrometer particles from the upper layers of the disk. 
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Table 1. Tensile strengths of selected materials. 
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are Coulomb unstable are denoted with a dashed line, whereas regions 
of the size distribution that are Coulomb stable are denoted by a solid 
line. 
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tion with atmospheric height. Regions of the size distribution function 
that are Coulomb unstable are denoted with a dashed line, whereas re¬ 
gions of the size distribution that are Coulomb stable are denoted by a 
solid line. 
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Fig. 8. Dust optical depth for an atmospheric plasma region where 
the Coulomb explosion of dust occurs. The top figure shows the ab¬ 
solute dust optical depth and the bottom figure shows the dust optical 
depth normalized by its maximum value so that a comparison of their 
profiles can be made. 
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